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We present an electronic circuit which simulates wave propagation in dispersive media. The circuit is an 
array of phase shifter composed of operational amplifiers and can be described with a discretized version of one- 
dimensional wave equation for envelopes. The group velocity can be changed both spatially and temporarily. 
It is used to emulate slow light or stopped light, which has been realized in a medium with electromagnetically 
induced transparency (EIT). The group-velocity control of optical pulses is expected to be a useful tool in the 
field of quantum information and communication. 



I. INTRODUCTION 

In an optical medium, the envelope of a light pulse travels 
at a group velocity v g , which is defined as the reciprocal slope 
of the dispersion relation for the medium. Recent advances in 
controlling the dispersive properties have enabled superlumi- 
nal or negative group velocities ll|,|2j], extremely slow group 
velocities |3], and even stopping light pulses QQl. 

In a classical textbook |6], Brillouin pointed out that the 
group velocity can exceed the light speed c in a vacuum 
and such a superluminal phenomenon does not conflict with 
causality. However, even after several experimental demon- 
strations fjl 

the superluminal light propagation has still 
invited some misunderstandings, because it is very counterin- 
tuitive that a light pulse runs faster than c. Recently we have 
proposed an electronic circuit which exhibits negative group 
delays. It produces an output pulse ahead of the input pulse 
dial- This circuit shares the same principle as light propaga- 
tion at superluminal or negative group velocities. The circuit 
is so simple that one can easily build it and examine how this 
puzzling phenomenon takes place. 

Extremely slow light has been realized by utilizing electro- 
magnetically induced transparency (EIT) |9]. In the EIT ex- 
periments, probe light and control light are tuned to a common 
excited level of A-type three-level atoms. The two light beams 
destructively interfere with each other through the common 
upper level and the light absorption disappears. The refractive 
index for the probe light, as well as the absorption coefficient, 
varies largely within very narrow bandwidth, called a "trans- 
parency window," where the absorption is vanished. Hence 
in the presence of the control light, the group velocity of the 
probe light becomes extremely small, compared with c. Hau 
et al. 1 3] succeeded in reducing the group velocity down to 
17m/s in laser-cooled atomic gas. They prepared a single- 
shot pulse with the duration of 2.5 /is and the pulse length of 
750 m in a vacuum, and fed it into the EIT medium. Owing to 
the slow down, the pulse length was compressed into 43 /im 
in the medium. 

The group velocity is proportional to the width of the trans- 
parency window, which is determined by the intensity of the 
control light; the weaker the control light, the smaller the 



group velocity. However, we cannot slow down the probe 
pulse unconditionally. The lower limit of the group velocity 
is given by the condition that the transparency window should 
be wider than the bandwidth of the probe pulse. Otherwise 
the pulse shape would be deformed by the absorption. This 
condition is often called the EIT condition. 

This limitation can be overcome by changing the intensity 
of the control li ght when the probe pulse is traveling in the 
EIT medium 1 10, 11]. The intensity of the control light is kept 
constant until the entire probe pulse fits into the EIT medium. 
Subsequently the control light is turned down smoothly, and 
the probe pulse is slowed down with its spatial profile un- 
changed. The pulse is stopped in the medium when the control 
light is eventually turned off. During the stopping process, the 
spectral width of the probe pulse is compressed and can be fit- 
ted in the transparency window, which narrows down accord- 
ing to the dimming of control light. The stopped pulse can 
be restarted by restoring the control light. The released pulse 
is identical with the original one, because all the processes, 
freezing and releasing, are carried out coherently. The char- 
acteristics of light pulses, such as the amplitude, phase, and 
beam profile, can be stored. Even the quantum state of light 
could be stored in the medium, and this method is promising 
as a quantum memory. 

This paper deals with the physics of the slow and stopped 
light, and introduces a circuit which simulates the light prop- 
agation with velocity control. This circuit model is very 
clear and instructive as that for the superluminal propagation. 
Moreover it is easy to construct the circuits and change the 
circuit parameters freely. We can set the time constant slow 
enough that the propagation can be seen as blinks of a series 
of light emitting diodes (LEDs), which monitor the output of 
each stage. The circuit behaves like a running message display 
of single row but the operation principle is quite different. 

In Sec. [II] we give a theoretical framework of pulse prop- 
agation in dispersive media and introduce the circuit for the 
simulation. In Sec.[ni]we show several experimental results, 
which correspond to slow light or stopped light experiments, 
and summary is given in Sec. lIVI 
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II. THEORY 
A. Pulse propagation and group velocity 

In a medium, the wave equation in one dimension with re- 
spect to the electric field, E(x, t), is given by 
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where P(x, t) is the polarization of the medium induced by 
the electric field. We introduce the envelopes, £(x,t) and 
V(x, t) as 

E(x, t) = £(x, t) e *(<"o*-faO + c.c, (2) 
P(x, t) = V(x, t) j(«>ot-kx) + c.c, (3) 

where c.c. stands for the complex conjugate of the preceding 
term. The carrier frequency luq is located near the center of 
the spectrum of the electric field. The associated wavenumber 
k will be determined later on. 

Assuming that £(x, t) and V(x, t) are slowly varying func- 
tions with respect to space and time, we can reduce the order 
of differentiations of each term of Eq. Q in the following 
way: 
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For slowly varying envelopes, the spectrum of the elec- 
tric field is confined within the narrow regions: [±a»o — 
5ui/2, ±wo + Soj/2], where Slu represents the spectral width. 
We have assumed that fractional changes of envelopes over a 
period or a wavelength are negligibly small. 

Then the susceptibility of the medium, x( w i x )> can be re- 
garded as a linear function of ui within the bandwidth and the 
relation between E (t, x) and V(t, x) can be written as 
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where xo = X^cOjXi = dx/du(uo)- The derivation of 
Eq. Q is given in Appendix |X] With the help of Eq. Q, 
Eq. (|6j can be represented as 



Considering the case of monochromatic fields (d£/dt — 
d£/dx = 0), we can determine the wavenumber k = nu) /c, 
where n = \/l + Xo is the refractive index. 

Then the former terms of Eq. (|9jl yield the wave equation 
for the envelope as 
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is the group velocity. 



B. EIT condition 

For EIT media, x( w ) changes sharply within the EIT win- 
dow, whose width is w. The slope of x is estimated as 
Xi = dx/du ~ a/w, where a is the variation of x in the 
window. Normally a corresponds to the optical thickness of 
the medium and can be of the order of unity. Assuming n ~ 1, 
we have an estimation 
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for the group velocity il 11 . Roughly speaking, w is equal to 
the Rabi frequency of the control light f} c , which is propor- 
tional to the amplitude of the control light 1 10]. The width 
w can be effectively reduced down to the atomic decoherence 
rate, which can be 100 kHz or even smaller. Therefore the ex- 
tremely slow velocity, v g ~ 10~ 9 c, or virtual stopping can be 
possible. 

The spectrum of the pulse must be contained within the EIT 
window. Namely, the spectral width Su> of the envelope must 
satisfy the condition: 



Slu < w. 



(13) 



This linkage between the group velocity and the bandwidth 
condition is a very important issue in slow light experiments. 
A similar restriction comes into play in cases of superluminal 
or negative group velocities QSl- 



C. Solutions 
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By substituting Eqs. ©, 0, and (|8j into the wave equation 
Q, we obtain 

o- j , w oXi\ d£ 2 d£ 
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If v g is constant in time and space, it is easy to show that 
the solution of Eq. fllOi is £(t — x/v g ). This means that the 
envelope propagates in the medium at the group velocity v g , 
keeping its shape. 

When v g is a function of x, the solution of Eq. dlOi can be 
written as 



£(t,x) = 
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where £(t, 0) = <fi(t) is the boundary condition at x — 0. The 
temporal profile of the pulse, or the spectrum, is conserved 
at any x, while the spatial profile shrinks (stretches) at the 
location of small (large) v g (x). This is the reason why the 
pulse length is compressed when the light pulse is fed from 
the vacuum into the EIT medium. 

Although the group velocity could be decreased further by 
reducing the width w of the transparency window, w must be 
wider than the bandwidth of the pulse, Sui, which is constant 
regardless of v g (x). Therefore, the spectral width Su> of the 
initial pulse sets the lower limit of the velocity. 

On the other hand, when v g is a function of t, the solution 
of Eq. i ll Ob is given by 



£ (t, x) = ip I x 



v g (t')dt' 
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where £(0, x) = i()(x) is the initial condition at t = 0. In this 
case, the spatial profile of the pulse is identical at any time t, 
but the temporal profile is changed according to v g (t). At time 
when v g (t) becomes small, the temporal profile, i.e. the pulse 
width, spreads out, and the frequency spectrum is compressed. 

If the EIT condition is satisfied initially, it is fulfilled all 
the time no matter how far we diminish the velocity v g or the 
transparency window w, because the spectral width Slu is de- 
creased accordingly. With this method, the stopping of light 
has been carried out. 



D. A circuit model 

We will show how to simulate the pulse propagation de- 
scribed by Eq. dlOi using an electronic circuit. First, we di- 
vide x into uniform steps as x = nAx (n is an integer). We 
discretize the second term of Eq. dlOt i as 
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where T = Ax/v g and v n (t) = £(nAx,t). We approximate 
the first term of Eq. dlOi as 



d£_ _^ d_ f v n+1 (t) +v n {t) 
dt ^ dt{ 2 

Then the discretized version of Eq. is obtained; 
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Introducing the Fourier transform of v n (t) 

V n {uj) = / V n (t) 



e~ luJt dt, 



we have the difference equation 

V n+ i(u) = H{u)V n {w) 



(19) 
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FIG. 1 : Elementary circuit for the simulation. 
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FIG. 2: Schematic diagram of overall setup and circuit parameters. 



where 



1 - iuT/2 
1 + iujT/2 ' 



(21) 



is the transfer function. 

The circuit illustrated in Fig.^h as the transfer function of 
Eq. d2 It . It is well known as an all-pass filter, which is used as 
a phase shifter 1121 . The circuit provides a group delay T for 
signals whose spectrum is limited within the bandwidth 1/T 
1 8]. Replacing the medium with the length Ax by the circuit 
with T = Ax/vg and cascading the circuits, we can simulate 
the envelope propagation described by Eq. dlOi . 

The quantity v& = 1/T, which appears in Eq. dl8> . corre- 
sponds to the velocity v g in Eq. dlOi . Even though the dimen- 
sion of v& is different from that of velocities, we hereafter call 
it "velocity." It is the number of stages passed by the pulse per 
unit time. 
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FIG. 3: Simulation of wave propagation at constant speed 
13/s. 
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TABLE I: Switch configurations and parameters. C = 82 nF. 



The arbitrary waveform generator (Tektronix: AFG320) 
provides Gaussian-like pulses to the input of the delay cir- 
cuits. The waveform is given by 



«o(t) = 



(t<0, i>5s), 

Voe -4(t-Uf/(5rf ( <i<5s), 



(24) 



where Vb = 1.0 V, t s = 2.5s, and St = 1.0s. The frequency 
spectrum of vo(t) is limited within Suo ~ 4/5r = 4.0 Hz. 



B. Constant speed propagation 

Figure[3]shows the experimental result for f<j = 13 /s. The 
Gaussian-shaped pulse propagates down at the speed of 13 
circuits per second with its shape unchanged. In this case, the 
spectrum condition, Eq. Q3> . is satisfied as 8lu/v& ~ 0.31 < 
1, which guarantees that the pulse travels with small distor- 
tion. The pulse length is estimated to be Sn = v&8t ~ 13. 



C. Pulse length compression 



E. Spectrum condition 

The above discretization is valid when the envelope £(x,t) 
is smooth enough that the changes of £(x, t) over Ax can be 
neglected. From Eq. (I18> . the fractional change between v n (t) 
and v n +i(t) is estimated as 
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where we have assumed v n (t) ~ v n+ i(t) and Su> is the band- 
width of v n (t). If the condition, 



Sui/vd = Stu ■ T <C 1, 



(23) 



is satisfied, the discretization is justified. This condition, Su> ■ 
T <C 1 is accidentally coincide with the condition that each 
all-pass filter works as an ideal group-delay circuit; H(uo) ~ 
exp(— iuiT). Thus we have found that the spectrum of the 
input signals should be restricted within 5ui <C v&. 



III. EXPERIMENT 

A. Circuit parameters and input signal 

The experimental setup is shown in Fig. [2] and the param- 
eters of the delay circuit are shown in the inset. Forty delay 
circuits are cascaded. Each of the output is monitored with an 
LED. The delay time T = 2RC, or the speed v& = 1/T, can 
be changed electronically with two analog switches (DG441) 
connected to the resistors. We show the relations between the 
switch configurations and the constants (R, T, and v&) in Ta- 
bled. 



In this section, we consider the situation where the speed 
v& is dependent on the location n. We divide the cascaded 
circuits into two regions; n < 25 (Region I) and n > 25 
(Region II). The velocity in each region is set as 

\v x (Region I; n < 25) 
I 1/2 (Region II; 25 < n) 

This case corresponds to the light propagation experiment 
where the pulse is fed from the vacuum (Region I) into the 
EIT medium (Region II). 

Figure [4] shows the result for the case of v\ — 13/s and 
V2 = 6.7/s(~ v\/2). At the boundary n = 25, the pulse 
propagation is decelerated to almost half the speed. It should 
be noted that the pulse has to be connected continuously at 
the boundary, which is marked by the dashed line, because 
the output of the 24th circuit is connected to the input of the 
25th circuit. Hence the pulse width St in time domain is kept 
unchanged for any location n, and the spectrum width Su> is 
also invariable. In the present case, for both of the regions, 
the spectrum condition is satisfied as Suj/v-i — 0.31 < 1, 
Slo/v2 — 0.60 < 1, and then the Gaussian shape is almost 
maintained through the propagation. On the other hand, the 
pulse length Sn in Region II is compressed into about half the 
length; Sn\ = v\St ~ 13 in Region I and Sri2 = v^St ~ 6.7 
in Region II. This is because the rear part of the pulse trav- 
els faster in Region I than the leading part runs in Region II 
after the leading part reaches the boundary. As mentioned be- 
fore, the same pulse compression occurs when the light pulse 
prepared in the vacuum is fed into the EIT medium. 

Next, we use the parameters: V\ = 13 /s, = 0.62 /s (~ 
fi/20). One might expect the pulse to be decelerated to 1/20 
and compressed spatially down to 1/20. But we obtain the 
result as shown in Fig. [5] In this case, the shape of the pulse is 
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(a) 3D 




FIG. 5: Experimental result for «i = 13/s in n < 25 and V2 = 
0.62/s in n > 25. The pulse is collapsed in n > 25, because the 
spectrum condition is violated. 
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(b) top view 

FIG. 4: Pulse length compression. Pulse which runs at vi — 13/s in 
7i < 25, slows down to V2 = 6.7/s in n > 25. 



collapsed completely in Region II. This is because the spec- 
trum condition is violated in Region II {5u}/v 2 = 6.4 > 1). 
The spectrum condition sets the lower speed limit: 8uj = 
4.0 /s. 

In the experiment of the slow light by means of EIT, similar 
spectrum condition exists. If the spectrum of the light pulse is 
wider than the transparency window of the EIT medium, the 
light is absorbed and the shape is strongly modified. 

D. Pulse freezing 

Here, we deal with cases where the speed v& is a function 
of time. For simplicity, we divide the time interval into three 
periods and in each period the velocity is kept constant; 

{v x (Period I; t < 4 s) 
v 2 (Period II; 4 s < t < 7 s) . (26) 
v 3 (Period III; 7 s < t) 

First we consider the case V\ — 13 /s and v 2 — = 6.7 /s. 
The result is shown in Fig. [6] At the end of Period I (t = 



4 s), the most of the pulse enters into the circuits. The pulse 
is then slowed down from the leading edge to the rear edge 
instantaneously. The pulse length Sn is unchanged, but the 
pulse width in time domain, St, doubles after the boundary. In 
other words, the frequency spectrum is compressed into half 
and Suj is halved. This means that the 6u> / v& remains constant 
(= 0.31). This situation corresponds to the case where we 
change the group velocity while the entire pulse is in the EIT 
medium (0 < n < 40). 

The violation of the spectrum condition, as discussed in the 
previous section, can be avoided by the strategy of changing 

in time. It is possible to attain much slower propagation or 
even to stop the pulse. We modify the time sequence as v\ = 
13 /s, v 2 — 0.62 /s, and vj, = 13/s. The result is shown in 
Fig-El ln Period II, the propagation of pulse is almost freezed, 
and released at t = 7 s. As mentioned above, the spectrum 
condition is always satisfied because Suo/v^ is conserved in 
any time. Hence, the pulse can propagate or be frozen without 
changing its shape in the time-varying system. 

In the EIT case, the change of the group velocity is at- 
tributed to the dark state polariton. The dark state polariton is 
a coupled mode of the electric field and the atomic coherence. 
The velocity of the polariton determined by the mixing ratio, 
which depends on the control light intensity. The atom-like 
polariton runs slower and the purely atomic polariton stays 
still. The pulse shape is imprinted in atomic states. The stor- 
age time is limited by the relaxation time of the atomic coher- 
ence. In the circuit model with a resistor connected parallel to 
each capacitor, we can simulate the decoherence. 

IV. SUMMARY 

In this paper, we introduced an electronic circuit which sim- 
ulates slow light and stopped light. Most of the properties in 
the simulation are common to that of the light propagation 
in the EIT medium. However the mechanism of the veloc- 
ity control is much more simple for the case of simulating 
circuit. Without knowing the details of EIT and dark state 
polaritons, we can understand the physics of slow light and 
freezed light. In addition, the simulation with the circuits has 
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FIG. 6: Changing the speed v g in time. Pulse which initially travels 
at vi = 13/s slows down to V2 = 6.7/s 




FIG. 7: Pulse freezing. 



the advantage that the parameters can be changed easily. This 
paper deals with three cases: constant v g , position dependent 
velocity v g (x), and time-varying velocity v g (t). It will be in- 
teresting to consider the general case where the velocity is a 
function of space and time, v g (x,t) (for example BOO . For 
cases where the implementation with EIT methods is difficult, 



the circuit model may likely be feasible. 

The present model deals with envelopes and can describe 
only one-way wave propagations. It cannot be used to repro- 
duce bidirectional phenomena such as reflections and interfer- 
ence, where both the forward and backward waves come into 
play. But it does reproduce one of the most fundamental prop- 
erties of waves; a disturbance propagates at a certain velocity 
keeping its shape through a medium. 
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APPENDIX A: POLARIZATION IN A LINEAR DISPERSIVE 
MEDIUM 

In this appendix, we derive the relation between the en- 
velopes £(t) and 'P(i). For brevity, we omit the explicit de- 
pendence on x. For linear media, the polarization P(t) in- 
duced by the electric field E(t) is given by 

/oo 
x(t-T)E(r)dr, (Al) 
-oo 

where x(t) is the response function. The dielectric suscepti- 
bility x(w) and the response function x(t) are connected with 
the Fourier transform: 

i r°° 

X(t) = ^J x(co)e^duj. (A2) 

The realness of x(t) implies x(— w ) = X*(w). Equation JA1I 
can be Fourier-transformed as 

P{lo) = e x{u)E(u). (A3) 

The Fourier transform of E{t) — £(t) e tuJat + c.c. can be ex- 
pressed as, 

/oo />oo 
£{t)e lUot e- lu)t dt + / £*{t)e- lUot e- lut dt 
-OO J —OO 

= £(uj - ojq) + S*(-w - w ), (A4) 
in terms of the Fourier transform of the envelope: 

/oo 
£{t)e- mt dt. (A5) 
-oo 

Assuming that the spectrum of the electric field is restricted 
around ±wo, we can use the Taylor expansion of the suscepti- 
bility xM: 

XM ~Xo+Xi(M -cj ), (A6) 
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where \o = x( w o)> Xi = dx/du>(u>o). When the losses in the 
medium are negligible, the susceptibility is real and must be 
symmetric; \(— w ) = x( w )- The moduli has been introduced 
in order to take into account the symmetry. 
From Eqs. (IA3I . (IA4> and (IA6> . we have 



P(w) ~ e xo^(^ - w ) + eoXi(^ - w )£ (w - cj ) 

+ eoxo^ - w ) + eoXi(-w ~ w )£*(-w 

and its inverse Fourier transform yields 

P(t) = eo fxo^(i) e Mot - *Xi^ e^ ot ^ + c.c. 



- UJ ) 

(A7) 



C.C. 



(A8) 



where we have used, 

£ «> - s 

Now we have Eq. 0: 



£(fi)e in *dfi, 
(zfi)£(fi)e if2 *dfi. 



7'(*)=eo(xo-txi^)f(t). 



(A9) 



(A10) 



It should be noted that only the second term contributes to the 
group velocity in deriving of Eq. dlOt . 
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